Demonstration of sufficient control for two rounds of quantum error correction in a 
solid state ensemble quantum information processor 
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We report the implementation of a 3-qubit quantum error correction code (QECC) on a quantum 
information processor realized by the magnetic resonance of Carbon nuclei in a single crystal of 
Malonic Acid. The code corrects for phase errors induced on the qubits due to imperfect decoupling 
of the magnetic environment represented by nearby spins, as well as unwanted evolution under 
the internal Hamiltonian. We also experimentally demonstrate sufficiently high fidelity control to 
implement two rounds of quantum error correction. This is a demonstration of state-of-the-art 
control in solid state nuclear magnetic resonance, a leading test-bed for the implementation of 
quantum algorithms. 



Introduction — One of the crucial requirements pQ for 
universal quantum information processing (QIP) is the 
ability to protect the fragile quantum information during 
computation - either by encoding the information in sub- 
spaces of the system's Hilbert space where it is protected 
from degradation by noise, or by an active scheme that 
detects and rectifies errors continuously or periodically. 
This latter, active, technique has been experimentally re- 
alized in liquid-state-NMR [2-4 and ion-trap [5 imple- 
mentations of a quantum information processor. In each 
of these cases, a quantum error-correction code (QECC) 
was used to protect against the particular errors present 
in the respective systems; and it was shown that even 
with imperfect encoding and recovery operations, em- 
ploying quantum error correction is advantageous. Obvi- 
ously, there is a limit to how many control errors can be 
tolerated before they overwhelm the error correction pro- 
tocol. Thus, the ability to demonstrate error correction 
is a highly relevant benchmark of high-fidelity coherent 
control. 

A natural question to ask is, whether one has high- 
enough-fidelity control to perform multiple rounds of er- 
ror correction, as would be required in a realistic compu- 
tation. Of course, to usefully perform multiple rounds, 
one needs a fresh supply of sufficiently pure ancillae to en- 
sure that entropy flows in the proper direction. However, 
assuming we have a fresh supply of ancillas, is it possible, 
with the current level of control, to perform meaningful 
multiple rounds of quantum error correction? 

In this manuscript, we report on the implementation of 
a three-qubit QECC that corrects phase errors induced 
by the environment, in a single-crystal solid-state NMR 
system. We also devise a way to experimentally deter- 
mine the entanglement fidelity of multiple back-to-back 
rounds of error correction, and use it to determine the 
entanglement fidelity of two rounds of the 3-bit phase 
code. In light of recent work on the characterization and 
control of such systems [SHE], as well as state initializa- 
tion [HI HO] , this current work signifies an advancement 



of one of the leading test-beds for QIP ideas. The rest of 
this manuscript is organized as follows. First, we describe 
the solid state NMR system, and the sources of noise af- 
fecting the qubits. We then describe the QECC imple- 
mented, and show the results from one and two rounds. 

System and error models — Building on the suc- 
cess of liquid-state NMR as a test bed of QIP ideas, 
Solid-state NMR systems offer [TT] intrinsically larger 
couplings, longer coherence times, the ability to pump 
entropy out of the system of interest into a spin bath 
and the potential for much higher initial polarizations. 
This comes at the cost of a more complicated Internal 
Hamiltonian, which makes the system harder to control 
in practice. 

The computational register under investigation is an 
ensemble of molecular nuclear spins in a macroscopic sin- 
gle crystal of Malonic Acid (C3H4O4). A small fraction 
(~ 3%) of the molecules are triply labeled with (spin- \ ) 
13 C to form an ensemble of processor molecules, spatially 
buffered from one another by molecules of the same com- 
pound but with natural abundance (~ 1%) Carbon nu- 
clei. During computation, the processors are decoupled 
from the 100% abundant spin- ^ protons in the crystal by 
applying a decoupling pulse sequence [H] to the protons. 
For this 3-qubit register, (dephasing) noise comes in the 
following forms: 

• Coherent phase errors due to pulse implemen- 
tation errors, phase transients, or unwanted evo- 
lution under the natural Hamiltonian (e.g. under 
the Zeeman term.) These are unitary errors (that 
cause no loss of coherence) and can therefore be in- 
verted if tracked properly, but in case tracking that 
evolution is not possible, quantum error correction 
becomes a valuable tool. 

• Incoherent phase errors due to Zeeman-shift 
dispersion or other inhomogeneities. The loss of co- 
herence is over the ensemble; each member of the 
ensemble sees a different value for some coupled 
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classical degree of freedom. Errors of this nature 
have been dealt with using refocussing techniques 
(e.g. spin echo), or by carefully designing the con- 
trol fields to generate the same evolution over the 
ensemble distribution of the inhomogeneous param- 
eter. Quantum error correction can be used with, 
or in lieu of, these other techniques to improve ro- 
bustness to ensemble errors. 

• Decoherent phase errors due to a coupling be- 
tween the system of interest and environment - 
an uncontrollable quantum degree of freedom- and 
loss of coherence occurs when this environment is 
traced over after the interaction. 

3-bit code — The 3-bit repetition code was introduced 
by Shor [13] as part of a 9-qubit code that is able to 
protect against an arbitrary single qubit error. The phase 
variant [TJ of the 3-bit repetition code encodes a single 
qubit in three qubits as follows: 
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where |±) = |0) ± |1), and the logical basis are denoted 
by {|0), |1)}. In the stabilizer formalism [El [16], the sta- 
bilizer group generators for this code are {XXI, I XX}. 
This code can be employed to correct for various sets of 
errors by choosing different decoding circuits - for this 
work, we design the decoding to correct for errors gen- 
erated by the set £ = {ZII , IZI , IIZ, III}. That is to 
say, with the same decoding circuit the code corrects a 
coherent selective phase rotation on one of the qubits, 

z e ._ e -i 6/2 zn = cos (0/ 2 ) HI- isin(0/2) ZII , 

and/or the dephasing map on one of the qubits, due to 
incohcrent/decoherent errors: 

A e (p) = cos 2 (0)lpl + sm 2 (6)ZpZ . 

A quantum circuit that accomplishes [14] the encod- 
ing, decoding, and error-correction steps is shown in 
Figure [TJb. The encoding process takes a qubit in the 
state a|0) + (3\1), as well as two ancillary qubits pre- 
pared in the |00) state, and outputs the 3-qubit encoded 

state a\+ + +) + (3\ ). After the error channel, the 

recovery process decodes the state on the information- 
carrying-qubit, and the other qubits carry syndrome in- 
formation about the errors that have occurred. The non- 
degeneracy of the code implies that each of the error ba- 
sis, in £, will leave a particular signature. It is straight- 
forward to show that syndromes 00, 10, 01, and 11 cor- 
respond to the occurrence of errors 77/, ZII, IZI, and 
IIZ, respectively. 

The figure of merit used herein to judge the perfor- 
mance of the code is the entanglement fidelity [17) . In 
particular, we use the expression for the single-qubit aver- 
age entanglement fidelity, which is experimentally acces- 
sible by measuring the fraction of surviving signal given 
input states X, Y, and Z [T8"| . 
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FIG. 1. Shown are the implemented quantum circuits for 
(a) Labeled-pseudopure-state (PPS) preparation procedure: 
a triple-quantum-coherence-filte (3QCF) is conjugated by a 
unitary operation that encodes (and decodes) the labeled 
pseudopure state j00X00|X in the triple quantum coherence 
|000><111| +|111)(000|; (b) The implemented quantum circuit 
of a three-qubit QECC, showing the encoding, decoding, and 
error-correction steps. The top two qubits are initialized to 
the 1 00) state, and the bottom qubit carries the information 
to be encoded. After the decoding and correction operations, 
the bottom qubit is restored to its initial state, while the top 
two qubits carry information about which error had occured; 
and (c) The procedure for two rounds: U p prepares X,Y, or Z 
inputs, and U s = {//, XI, IX, XX} toggles between the dif- 
ferent syndrome subspaces. i.e. The experiment is repeated 
four times, cycling through the different U s , and then the 
results are added, similar to a standard phase cycling proce- 
dure. 



Experiment — The experiments were performed 
in a static field of 7. IT using a purpose-built NMR 
probe. Shown in Figure [2] is a proton-decoupled 13 C 
spectrum, following polarization-transfer from the abun- 
dant protons, for the particular orientation of the crys- 
tal used in this experiment. A precise spectral fit gives 
the Hamiltonian parameters (listed in the inset table 
in Figure [2]), as well as the free-induction dephasing 
times, TJ, for the various transitions; these average at 
~ 2ms. The dominant contribution [7] to T| is Zeeman- 
shift dispersion, which is largely refocused by the control 
pulses. Other contributions are from intermolecular 13 C- 
13 C dipolar coupling and, particularly for C m , residual 
interaction with neighboring protons due to imperfect de- 
coupling. The carbon control pulses are numerically opti- 
mized to implement the required unitary gates using the 
GRAPE [H] algorithm. A typical pulse is 1ms long, and 
is designed [20] to have an average Hilbert-Schmidt fi- 
delity of 99.8% over appropriate distributions of Zeeman- 
shift dispersion and control-fields inhomogeneity. 

The system is initially prepared in the labelled pseu- 
dopure state (PPS) [2iTf2"3] -expressed in the product op- 
erator formalism- Pi — \ [HI + e{I + Z)(I+ Z)X] , where 
e ~ 10~ 5 . The completely mixed component is ignored 
for the rest of the discussion, as it does not participate in 
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FIG. 2. Malonic acid (C3H4O4) molecule and Hamiltonian 
parameters (all values in kHz). Elements along the diagonal 
represent chemical shifts, uji, with respect to the transmitter 
frequency (with the Hamiltonian ^2 i ntOiZi). Above the di- 
agonal are dipolar coupling constants (X/i<j ^Dijip, ZiZj — 
XiXj — YiYj), and below the diagonal are J coupling con- 
stants, (X^i<j ^Ji,j(ZiZj + XiXj + YiYj). An accurate nat- 
ural Hamiltonian is necessary for high fidelity control and is 
obtained from precise spectral fitting of (also shown) a proton- 
decoupled 13 C spectrum following polarization-transfer from 
the abundant protons. The central peak in each quintuplet is 
due to natural abundance 13 C nuclei present in the crystal at 
~ 1%. (for more details sec [7, 10 and references therein.) 



the unital dynamics. This preparation is achieved by con- 
trol pulses and phase cycling (temporal averaging), and 
can be thought of as a projection along (I + Z)(I + Z)X . 

As shown in Figure [l}a, the phase cycling actuates a 
triple- quantum- coherence filter (3QCF) |24| by exploit- 
ing the n-proportional phase acquisition of n-coherence 
quantum states under Z-rotation. And conjugating the 
3QCF with transformations that encode (and later de- 
code) the (/ + Z)(I + Z)X coherence in the triple- 
quantum-coherence, 1 000^1 1 1 1 + |111|000|, realizes an 
effective projector on the labelled-PPS. 

We first examine the performance of the 3-bit phase 
code under the natural evolution of the system: between 
the encoding and recovery operations, the system is left 
to evolve unobstructed under the full natural Hamilto- 
nian, both homonuclear and heteronuclear parts. Exag- 
gerated as it is, this is a useful test of the code's abil- 
ity to correct for coherent errors from uncertainties in 
the natural Hamiltonian, or imperfect decoupling of the 
magnetic environment. The experimentally determined 
entanglement fidelities are shown in Figure [3j demon- 
strating the advantage of quantum error correction. The 
syndrome information (inset in Figure [3| indicate that 
the dominant phase error is on the methylene carbon, 
C m . The non-montonicity of the unencoded and de- 
coded data indicate that the error is, at least partially, 
coherent. However, full simulation of the dynamics of 
the carbon subsystem suggest a longer timescale for the 
homonuclear coherent effects. Moreover, the timescale 
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FIG. 3. Experimentally determined entanglement fidelities 
for unencoded (blue squares), before (green diamonds), and 
after (red circles) the correction step of one round of QEC. 
After encoding, a variable delay is implemented before the 
recovery process. During the delay, the system evolves un- 
der the heteronuclear and homonuclear terms in the natural 
Hamiltonian. The non-montonicity in the unencoded and de- 
coded data are indicative of the presence of a coherent error. 
The inset shows the intensities measured for the different syn- 
dromes; the dominant error is a phase rotation on the bottom 
qubit (C m ). 



of the revival of the signal is consistent with the cou- 
pling strength between C m and the methylene protons, 
which leads us to conclude that this coupling is respon- 
sible for the non-motonicity in the entanglement-fidelity 
decay. This conclusion is further supported by the fol- 
lowing results, where this coupling is partially averaged 
using a heteronuclear-dipolar decoupling pulse sequence. 

Next, the 3-bit phase code is employed to protect a sin- 
gle qubit against errors from evolution under the natural 
Hamiltonian of the carbon subsystem as well as residual 
heteronuclear couplings between the carbons and pro- 
tons due to partial decoupling of the protons using the 
SPINAL64 sequence [H] at amplitude of 70kHz. From 
the syndrome information shown in Figure |4j the major 
contributions are from phase rotations on Ci and C m . 
This is to be expected, since, in this orientation, and in 
this rotating frame, the Zeeman shifts of these two spins 
are the dominant terms in the internal Hamiltonian. 

Two rounds — We devise a way to experimentally 
determine the entanglement fidelity of multiple rounds of 
error correction and use it to experimentally determine 
the entanglement fidelity of two rounds of the 3-bit phase 
code. After the first round of error correction, the sur- 
viving polarization from the various input states is dis- 
tributed over the various subspaces of the Hilbert space 
corresponding to the various syndromes. For the second 
round, for each syndrome, we project into the subspace of 
the syndrome and perform error correction in that sub- 
space, and then sum over all possible syndromes. For 
each syndrome, this projection is implemented as a trans- 
formation (denoted by U s in Figure [lie) that swaps the 
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rounds of correction performed. In this sense, the appli- 
cability of the scheme is very limited, but it is sufficient 
for our purposes. 

Conclusion — We were able to demonstrate the ad- 
vantage of performing quantum error correction to pro- 
tect against relevant, naturally occurring phase errors 
-coherent, incoherent and decoherent- that arise in a 
solid-state system. We have shown that this is possible 
by achieving state-of-the-art control on a 3-qubit system. 
Moreover, we have shown that with these control fideli- 
ties, multiple rounds of QEC are possible. This is partic- 
ularly significant in a system where it has been previously 
shown that entropy can be efficiently removed from the 
system of interest to the environment [HI [TU] . 
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FIG. 4. Summary of experimental results for the partial de- 
coupling map: the system evolves under the natural Hamilto- 
nian as well as 70kHz decoupling fields that partially modu- 
late the heteronuclear interactions (between the carbons and 
protons.) Shown (in top figure) are the single-qubit entan- 
glement fidelities in the cases where no encoding is employed 
(blue dots); or one round of the 3-bit code (red crosses); or 
two rounds of the 3-bit code (black asterisks), where the inter- 
action interval is split to two equal intervals. The dashed lines 
are quadratic fits to the data, and are included to guide the 
eye. Also shown (bottom) is the signal after one round of er- 
ror correction as distributed over the various error-syndrome 
subspaces. In this case, the dominant errors are phase flips 
on the top and bottom qubits, which are encoded on Ci and 
Cm, respectively. 



information in that subspace with the subspace where the 
ancillas are in (I+Z)(I+Z) -or |00)(00|- and then project- 
ing unto the latter subspace using the same protocol for 
initial state preparation. The quantum circuit describing 
the protocol is shown in Figure[T]-c, and the experimental 
results are shown in Figure [I] 

The results show that, for long interaction intervals, 
there is an advantage to performing two rounds of error 
correction with our current level of control. The initial 
drop (at zero interaction interval) in the experimentally 
determined two-round entanglement fidelity is mainly 
due to the projection operation, which is not needed if 
pure ancillaeare available. 

The scheme requires a number of experiments that 
grows as s m_1 , where s is the number of possible nondc- 
generate syndromes of the code, and m is the number of 
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